Solution of the fermionic entanglement problem with interface defects 

Viktor Eisler^ and Ingo PescheP 

^ Niels Bohr Institute, University of Copenhagen, 

Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

"^ Fachhereich Physik, Freie Universitdt Berlin, 

Arnimallee 14, D-14-195 Berlin, Cermany 

We study the ground-state entanglement of two halves of a critical transverse Ising chain, 
^^ ' separated by an interface defect. From the relation to a two-dimensional Ising model with 

a defect line we obtain an exact expression for the continuously varying effective central 
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charge. The result is relevant also for other fermionic chains. 



INTRODUCTION 



In one-dimensional systems, already a local perturbation can have important and interesting 
effects. This holds for the transport properties, but also for the entanglement between two pieces 



of a quantum chain. The influence of a defect located at the interface between them was first 

n -■ 



<— > 

j^ ! investigated by Levine for a Luttinger liquid [1]. In this case, also studied on a lattice in [^, the 

interaction leads to a strong renormalization of the defect strength |3|] and thus to special properties. 
^ , The situation is simpler for a fermionic hopping model, or XX chain, which was investigated in [J]. 



In this case it was found that the logarithmic dependence of the entanglement entropy S on the 
subsystem size L remains intact, but the prefactor depends continuously on the defect strength. 

cn 

K*" ' Thus, for I' interfaces 

■^ 

^; S = u^lnL + k (1) 

in 

f^ ' with a variable coefficient which replaces the central charge c of the pure system. The dependence 



of Ccff on the defect was determined rather precisely by numerical calculations, both in XX chains 
and in transverse Ising (TI) chains [5[. But although a simple analytical expression fits the data 
qualitatively, no exact formula has been given so far. 

In this paper, we want to look at this problem once more and present an analytical treatment 
based on the connection with a two-dimensional classical system with a defect line. For the TI 
quantum chain, this is a two-dimensional Ising model where defect lines are in fact marginal 
perturbations which lead to variable local magnetic exponents and have been studied over the 



years by various methods, see 



6| for a review and [7\ 



for a review and [7] for a more recent approach using boundary 



conformal field theory. Our techniques will be simpler and purely fermionic. Using conformal 
mapping, we will reduce the considerations to the study of the transfer matrix in a strip with 



defect lines parallel to the edges 



B. 



From its single-particle excitations we will obtain not only the 
spectrum of the reduced density matrix, but also the formula for the effective central charge. Our 
approach is similar to a recent calculation by Sakai and Satoh, who treated conformal interfaces 
between two different bosonic systems with c = 1 [9]. Actually, we were motivated by this work, 
because the formula found there does not describe the results found in the fermionic XX and TI 
chains. Nevertheless, our final result will turn out to be closely related. 

In the following Section 2 we set up the problem and present some numerical results for the 
spectrum of the reduced density matrix for later comparison. In Section 3 we describe the connec- 
tion of the TI chain with the two-dimensional Ising model, including the necessary renormalization 
of parameters. In Section 4 we determine the proper transfer matrix in the strip and its single- 
particle spectrum. This allows us to find the entanglement entropy in Section 5 and to give a closed 
formula for Cefj . En route, we also obtain the largest eigenvalue of the reduced density matrix. In 
Section 6 we sum up our findings and discuss the extension to other chains and defects. Additional 
details of the transfer matrix calculation are given in Appendix A, while in Appendix B we show 
that with our spectra one can also rederive the local magnetic exponent, which is a check on the 
approach and a further application. 

II. CHAIN PROBLEM 

In the following we will study the transverse Ising chain with 2L sites, open ends and Hamilto- 
nian 

n=-L+l n=-L+l 

where the o"^ are Pauli matrices at site n. We are interested in the entanglement of one of the 
half-chains with the other in the ground state. To have a modified bond at the interface, the 
parameters are chosen as hn = h, Jn = J for n ^ and Jq = t J. Similarly, a modified transverse 
field at the interface is described hy hn = h for n ^ and ho = th. We work at the critical point 
given hy J = h = 1. Then H is normalized such that the velocity of the single-particle excitations 
in the homogeneous system is v = 1 and only one parameter, the defect strength t, remains. 

To calculate the entanglement, one introduces fermions in ([2]) via the Jordan- Wigner transfor- 
mation and then determines the fermionic single-particle eigenvalues ei of the operator 71 in the 



exponent of the reduced density matrix 



P = ^eM-n) (3) 



from the correlations [10]. The entanglement entropy S = — Tr(pln/9) is then given by 

^ = EMl + e--) + 5:^ (4) 

I I 

Numerical calculations along these lines were done in [5| for the case of a ring with two equal 

defects on opposite sides. This configuration has essentially the same features as the open chain 

and one obtains the result ([TJ with v = 2. The effective central charge Ceflf(t) rises from zero for 

t = (where the system is cut in two) to one-half for f = 1 (where it becomes homogeneous). It is 

shown in Fig. 2 of Ref . [5| . 

The same function had been found before, also numerically, for the case of a segment in an 

infinite XX chain with a hopping defect at one interface [^. The other, unmodified interface 

then contributes an additive constant to Cefj (t). The background of this universality is the close 

connection between the two kinds of chains. By taking a critical TI chain with a bond defect and 

another one with a field defect, both of strength i, superimposing them and using dual variables, 



one arrives at 
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1^ 



H = -\ E(«+i + ^K+i) - \ * (^oVf + olaX) (5) 

n 

It was shown in [13] that the XX system with 2L sites in the subsystem and the combined TI 
systems with L sites have the same set of eigenvalues £;. This holds even in the inhomogeneous 
case and can easily be checked numerically. As a consequence, SxxC^L) = S^j{L) + S^j{L), which 
explains the common Cgfj. 

Given the central role of the single-particle spectrum oi Ti, it is important to see how this 
spectrum varies with the defect strength. This is shown in Fig. [1] for TI chains with 2L = 300 
sites. The characteristic feature is the development of a gap at the lower end of the spectrum and 
an upward shift of the whole dispersion curve, which increases roughly like ln(l/t) as t goes to zero. 
This upward shift causes a decrease of S for a fixed value of L. In particular, it makes S vanish if 
t = and the system is cut in two. The logarithmic variation of S with L, on the other hand, is 
related to a gradual lowering of the eigenvalues as L increases. These features will be found again 
in the analytical treatment, to which we now turn. 
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FIG. 1: Single-particle eigenvalues ei as a function of the defect strength for TI chains with 2L = 300 sites. 



III. RELATION TO A TWO-DIMENSIONAL PROBLEM 



In the case of a homogeneous chain, the TI Hamiltonian commutes with the diagonal transfer 
matrix of an isotropic two-dimensional Ising model and its ground state can be obtained from the 

id . Il4l | . The reduced density matrix is then given by a strip 



partition function of a long Ising strip 



with a perpendicular cut as shown in Fig. [2] on the left. One expects that the chain with a defect 
is related to a planar Ising system with a defect line as indicated. However, the relations in 15| 
show that by diagonally layering the Ising lattice one cannot obtain a commuting TI Hamiltonian 
with only a single modification (bond or field). Therefore we consider a lattice in the normal 
orientation and its row or column transfer matrices. It is well-known, that in a proper anisotropic 
(Hamiltonian) limit, they are just the exponential of a TI operator. This also holds with defect 
lines. For the case of a ladder defect, we can refer to Fig. [3l If the horizontal couplings are large 
and the vertical ones are small, the column transfer matrix, running in the direction of the ladder, 
is given by 



T = A exp{-2K^ H) 



(6) 



with // as in ([2]) and A a constant. The quantity K2 is the dual coupling of K2, sh(2i^2) = 
l/sh(2i^2)i and thus also small by assumption. The parameters in H are /i„ = 1, J„ = K1/K2, 
Jq = Kq/Ki Jn- At the critical point K2 = Ki and one obtains the critical TI model with a bond 
defect given by t = Kq/Ki. 

This is, however, not quite enough, since we want to use a conformal mapping in the following, 
for which the 2D system should be isotropic. In the homogeneous case, one can achieve this by a 




FIG. 2: Representation of p for a half-chain by two-dimensional partition functions. Left: Original repre- 
sentation. Centre: Simplified annular geometry. Right: Strip geometry obtained by mapping. The defect 
line is shown dotted. 



rescaling such that the velocity of the excitations becomes v = \. Thus one omits the factor 2K2 
in ([6]) and has T = Aexp(— f/^). This is then valid in the continuum limit or for large distances. 
With a defect, the situation is more subtle. The local magnetic exponent /3i, which also gives the 
power-law decay of the correlations parallel to the defect, is for the ladder [y, [la, ll7| 



(5i = —^ arctan 



vr^ 



thi^i 



(7) 



thi^o, 

where the couplings are taken at criticality. Thus it matters whether one is in the Hamiltonian 
limit, where the argument in the arctan is Ki/Kq, or in the isotropic case, where it is thi^i/thivTo 
with thKi = th i^c = V^ — 1- Therefore the TI Hamiltonian only corresponds to an isotropic 
system if also the defect parameter t is renormalized to 

thKo 



t 



bond defect 



(8) 



thi^i' 

For a chain defect with modified horizontal couplings Kq, analogous considerations apply. The TI 
Hamiltonian then contains a modified field given by Kq/K2 where now the dual couplings enter. 
These also appear in the local exponent and one has to renormalize the defect parameter to 

thET, 



thi^*' 



field defect 



(9) 



where one can use K2 = Ki = Kc- 

Having related the TI chain to an isotropic Ising lattice, we can now turn to the partition 
function. Instead of studying the geometry shown on the left in Fig. [21 we will consider the one in 
the centre, where the system is an annulus with a small inner radius a of the order of the lattice 
constant and outer radius R = La. This could be achieved approximately by the mapping used in 



18l |. The pieces in the lower and upper half-plane can now be viewed as azimuthal corner transfer 
matrices (CTMs) spanning an angle of 180° and containing a defect line in their centre. Using the 
mapping w = \nz, the annulus becomes the strip shown on the right of Fig. [2] with width \n{R/a) 
and height 2tt. One can now calculate the partition function Z{n) for a strip formed by repeating 
n such units to give a height of 27rn, and closed in the vertical direction. Then S follows from the 



usual formula, used also in ^^ 

S= (l -^]ln Z(n) . (10) 

However, since the transfer matrix needed to calculate Z{n) is just the transform of the corner 
transfer matrix which upon squaring gives /9, one can also obtain the single-particle eigenvalues in 
% directly from it and then use (H]). 

IV. TRANSFER MATRIX 

We now consider the strip and discretize the problem again by inserting a lattice with 2M rows 
and A^ columns. The couplings are isotropic, but we leave them anisotropic for the moment. The 
basic unit is half the strip with M rows and one defect line in the middle. This is shown, for the 
case of a ladder defect, in Fig. [3l The row transfer matrix for the whole unit will be called W . It 
has the form W = V^'^''^VqiV^'^''^ where V = V^ V^V^ is the symmetrized row transfer matrix 
for one step and Vqi describes the ladder. The quantities V\ and V2 contain the vertical and the 
horizontal bonds, respectively, and are given by 

(11) 

where we left out a prefactor in Vi which is not important here. The second, fermionic forms arise 
after a spin rotation and a Jordan- Wigner transformation. The ladder operator is Vqi = VqV^i^ 
where Vq has the same structure as Vi with Kq replacing Kf. 

Although such product matrices as V and W can be diagonalized for the open boundaries we 
want to treat [ly], the procedure is much simpler for periodic boundary conditions. We therefore 
assume these for the moment, which will give us the functional form of the eigenvalues. A Fourier 
transformation Cn = ^~^'^ ~7^ Yliq^^'^^^q then gives 

Vl'^ = n exp [kI{c\c, + ct_,c_,)} 

(12) 
^2 = n ^^P Y^'i cos q{c\cq + cl^c-g) - smq{clc^_g + c^qCg) | 

q>0 
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FIG. 3: Lattice and couplings in the strip geometry. 

where the momenta q are different for even or odd total fermion number. The problem now 
separates in (g, —q) subspaces, and is solved by introducing new operators aq,l3q via a Bogoljubov 



transformation 



2C| 



a„ 



cos (fq — Sm Lfq 

sin ipq COS ipq j \ 0q 

The diagonal form of the operator V then reads, after changing to k = tt — q 



y = Aexp ( X]7fc(al«fc + 4/3fc) I 
\fc>o / 



with the well-known dispersion of the homogeneous system 



(13) 



(14) 



ch 7fc = ch 2KI ch 2K2 - sh 2Kl sh 2K2 cos k 



(15) 



In a second step, one uses this result to diagonalize W . Some details are given in Appendix A. W 
has the same exponential form ()14p as V but with other operators and single-particle eigenvalues 
uJk given by 



ch Uk = ch 2 A ch M^^ + sh 2 A sh M^j^ cos 2^9^ 



(16) 



where A = K^ — K^. This is completely general. We now specialize to an isotropic critical system 
and small k. Then one has 7^ = k, cos 2ipk = k/^/2 and the second term in ()16p can be neglected. 
If one writes Mk = e^ for the eigenvalues of the homogeneous system, the relation becomes 



ch LVk = ch 2 A ch e^ 



(17) 




FIG. 4: Dispersion relation for the single-particle excitations in W for several values of the defect strength. 



The basic feature is that the defect introduces a gap 2A into the spectrum, because it makes the 
system locally non-critical. The parameter ch 2A can be rewritten as 



eh2A = ir^ + "'^' 



2 \thKi thKo, 
which gives via ([8]) a very simple relation with the chain parameter t 



, . 1 / 1 
ch 2A = - t + - 

2 V t 



(18) 



(19) 



In Fig. S]we show the resulting dispersion curves for several values of t. One sees that they become 
rapidly linear for larger values of e and show an upward shift, which from pzp is Inch 2 A and 
thus varies logarithmically with t for small t. These are exactly the features found in Fig. [1] for 
the numerical eigenvalues of the operator Ti. According to the remark at the end of section 3, 
these should correspond to 2uji., since one needs two transfer matrices to build the full system. A 
comparison of the numerical values of 2u}{t) and 2oj(t = 1) verifies (|17p relatively well, but there 
are shifts which may have to do with the open boundaries. The allowed values of the momenta k 
in this case have to be determined from an equation which is known in the homogeneous case [19l | 
and will be modified by the defect. In particular, they will vary with the defect strength. 

So far our consideration has been for a ladder defect in the lattice. The case of a chain defect 
can be obtained very simply by making a dual transformation in the transfer matrices Vi and V2. 
This leaves the spectrum invariant, but changes vertical bonds into horizontal ones. The resulting 
chain defect then has coupling Kq. Changing this into the desired Kq, one finds that ch2A has 
again the form ()19p with t now given by ([9]). Thus bond and field defects in the TI chain lead to 
the same transfer matrix spectra. 



One should mention that a similar situation was found for CTM spectra in critical Ising models 
with a radial perturbation of the couplings decaying as 1/r [6|, |21]. This problem can be mapped 



to a homogeneous strip which is slightly non-critical. The dispersion relation then assumes the 



relativistic form w^ = ^/A^ + e|. Our result reduces to this form for small A. 
For the following considerations we define the parameter 

sin(2 arctan t) = sin { 2 arctan — ] (20) 



ch2A ' 'V t 



which will turn out to be the quantity used in Ref. 



V. ENTANGLEMENT ENTROPY 

We are now in the position to evaluate S. For this we insert 2^;^ as single-particle eigenvalues 
in the expression (JH) and assume a large value of A^, so that the sum can be converted into an 
integral via 

^ _ / dk = 4i- d£ (21) 

The gradual shift of the allowed values of k for free boundaries does not play a role in this limit. 
The mapping gives 2M/N = 27r/lnL, so the prefactor of the second integral is InL/vr^ and we 
have the formula 

S=^nnL (22) 

with 

I = h + h= deln(l + e-2-)+ / de ^— — (23) 

Thus the entanglement entropy varies logarithmically and the coefficient is given by the integrals 
in /. Inserting chw = che/s, it depends on the single parameter s which in turn is determined by 
the defect strength t. In principle, this is the complete solution of the problem. 

In the homogeneous case where t = s = 1, the two integrals both have the value 7r^/24 and one 
obtains the standard result S = c/61nL with c = 1/2. For general values of s, the integrals cannot 
be evaluated, but only brought into a simpler form. Differentiating twice with respect to s and 
substituting x = ch e/ s gives 



1 r 



dx 



/ 



l/s wx"- - s 



arch X 



\/x^ — 1 ^ 



(24) 
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FIG. 5: Effective central charge Ccff(i) as a function of the defect strength t and the normaUzed quantity 

I21is)/Tr\ 



This can be evaluated using partial integrations with the simple result 



1 



Ins 



I"is) 



Integrating this again with the boundary values 1(0) = /(O)' = then leads to 



(25) 



I{s) 



1 + s) ln(l + s) + (1 - s) ln(l -s) Ins + (1 + s)Li2(-s) + (1 - s)Li2(s) } (26) 



where Li2(s) denotes the dilogarithm function defined by 

ln(l — x) 



U2{z) 



dx- 



(27) 



/o X 

The function 12/7r^I(s) is shown in Fig. [5] and rises smoothly from to 1 as s increases. At s = 0, 
it is non-analytic and varies as (s^/2) ln(l/s). The quantity Ces = G/vr^/ as a function of the defect 
strength t is also shown in Fig. [S] and varies between and 1/2. It has the same non-analyticity 
near t = 0, found already in [5], but varies quadratically near t = 1. If one takes the numerical 
data of [51], they lie perfectly on the curve. The same holds for the data on the XX model. Since 
s{t) = s{l/t), one obtains the same c^s for weakened or strengthened defect bonds. This was also 
observed in the numerics. According to the previous findings, field defects will also lead to the 
same results. 

A quantity closely related to S is the largest eigenvalue wi of the reduced density matrix p. It 



plays a role in the so-called single-copy entanglement problem 
quantity Si is given by the first term in / 

5*1 = — ^-/i In L = —:— In L 



22t |. If one writes Si = —Inwi, the 
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(28) 
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FIG. 6: Coefficient of the single-copy entanglement Kcs{t) as a function of the defect strength t and the 
normalized quantity 12/i(s)/7r^. 



with a coefficient k^s analogous to Ces- Therefore wi varies as a power of L and the exponent 
is — Keff/6. In the homogeneous system, k = c/2 because then the integrals /i and I9, have the 
same value. This is a general result for (homogeneous) conformally invariant systems |23l425|. For 
arbitrary s, one can calculate /i in the same way as /. For the second derivative one finds 



I'lis) 



1 



1 



+ 



1 

2^2 



ln(l 



(29) 



which upon integration gives the simple result 



1 



/i(s) = - Li2(s)+Li2(-s) 



(30) 



The quantity 12/tt'^Ii{s) is shown in Fig. E]and seen to rise from zero to 1/2 with a vertical slope at 
s = 1. In Keff (t), the slope is zero, but the function is non-analytic, varying as (1 — t)^ ln(l/(l — t)) 
near t = 1. For small t, the behaviour is Gi^/vr^, which was found already in [5] via perturbation 
theory. In [5|] the whole function Kcff(0 was also calculated numerically and again the data match 
the formula given here. 

Finally, we want to make contact with the result found by Sakai and Satoh in J9l]. There, in 
Eq. (4.6), a c^q was given which consists of a term linear in a parameter s minus an integral of 
bosonic origin which seems very different from the fermionic I{s). However, its second derivative 
is exactly ()25p and it actually equals /(s) as one can also verify by numerical integration. But it 
enters with a negative sign, and the additional linear term gives a finite slope at s = 0. Therefore 
the given expression cannot apply to our situation. 



12 



VI. DISCUSSION 

Our calculations have given the asymptotic behaviour of the entanglement entropy between two 
halves of TI and XX chains separated by interface defects. In particular, the effective central charge 
was found directly in terms of the corresponding defect parameter t. The same holds for Kefj related 
to the largest density-matrix eigenvalue. The continuous variation of both quantities is analogous 
to that of the magnetic exponent /3i, and this analogy becomes even closer if one determines /3i as 
in Appendix B. Thus all of them reflect the marginal character of the perturbation. 

Physically, one expects that these features are correlated with more elementary properties of 
the problem. For the XX chain, the obvious one is the transmission through the defect. This 
can be expressed through the two scattering phases one obtains from the S matrix. At the Fermi 
surface for half filling (kp = '7r/2), these are given by 



B 



5± = ± (2arctani- -") (31) 



and were already used in [J] to guess an approximate expression for Ceg . The transmission coefficient 
then is 



2 i'5+ - 6- 



T = cos^ ( -!— j = sin^(2 arctan t) = s^ (32) 

which gives a direct physical interpretation to the parameter s. The TI model is not as straight- 
forward, but one finds the same result. Such formulae also appear in the treatment of conformal 
defects [26[]. They open the possibility to apply the results found here also to other situations. 
The simplest case is the XX chain with arbitrary filling. Then it can be checked numerically that 
the corresponding Tikp) gives the correct Ces- But one can also consider more complicated types 
of defects, where a direct calculation as done here would be quite difficult [22]- An interesting 
problem would also be the extension to a quench from a homogeneous system to one with a defect. 



The numerics show, that Ces then appears in the time dependence 



B. 
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APPENDIX A. 

In this appendix, we give some more details on the diagonahzation of W for completeness. Such 
calculations have, of course, been done also in other studies dealing with layered Ising lattices, see 



e.g. 



2mm 






(33) 



We begin with V. Forming Heisenberg operators with Vi and V2, one has 

1/2 ^ ( e-^i 
\ e-^'i 

C2 — 5*2 cos q S2 sin q 
5*2 sin q C2 + S2 cos q 

with the notation Cj = ch.2Ki and Si = sh2i^j for i = 1,2. This gives for V 

_ / e-2^i (C72 - S-a cos g) S'asing j I ^^ 1 (u) 

\ S2smq e^^t{C2 + S2Cosq) ) 

The matrix on the right has determinant 1 and thus eigenvalues e '^'' where 

ch 7g = CI C2 + SI S2 cos q (35) 

It is diagonalized by the canonical transformation (J13p where the angle is, for the isotropic system 
at the critical point 




2cos2f ^ ^ 

COS 2^, = J ^^^^ (36) 



The matrix V then has the diagonal form given in ()14p . 

In the second step one needs the Heisenberg operators of aq and jSq with V '"^ and Vqi. These 
are 




"9 I y-M/2 = I 

y e^^«/2 

/ c2e-2 

^01 1 : I ^ol' 




_^ , c2e-2A + sV^ cs(e2^-e-2^) 




(37) 



where A = i^Q — K\ , c = cos 9?^ and s = sin Lpq . This finally gives for the transfer matrix W 
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with 

a = (ch 2A - sh 2A cos 2(pq)e-^^^'' 

b = (ch 2A + sh 2A cos 2(/jg)e^^« (39) 

d = sh2Asin2(^g. 

Again the determinant equals 1 and the eigenvalues have the form e «. This leads to the result 
(fT6]) in the text. 

APPENDIX B. 

The local magnetic exponent /3i was used in section 3 to find the correct corespondence between 
the defect parameters in the chain and in the Ising lattice. If /3i were not known, one could also 
have calculated it with our approach. In the CTM method, order parameters are obtained by 
fixing the spins at the outer boundary and then calculating the expectation value niQ of the central 
spin. For free fermionic systems, this gives rriQ as a product involving the single-particle CTM 
eigenvalues [6|]. Taking the logarithm and converting the sum into an integral via ()2ip one finds in 
our notation 

1 f°° 
In mo = In L — / de In th uj (40) 

^ Jo 
The coefficient of InL is then — /3i. The same formula was used also in [2u\. In contrast to I{s), 
this integral can be evaluated in closed form. Substituting z = 1/che one has 

13^ = -— I ^^— -^ [ln(l + sz) + ln(l - sz)] = ^ arcsin^ s (41) 

Since s = sin(2 arctan t) = sin(2arctan j), one has to make a choice here. Taking the second form 
gives 

/3i = — arctan^ i-\ (42) 

which is the ladder result ([7]) if one writes it in terms of the lattice parameters. However, since 
(|4ip cannot exceed the value /3i = 1/8 of the homogeneous system, this expression is valid only for 
t > 1, or Kq > Ki. For f < 1, one has to take the first form for s. But in this case, when the system 
is slightly in the disordered phase, the fixed boundaries lead to an additional low-lying eigenvalue 
ujQ wich vanishes as a power of L and thus contributes to /3i . Such a situation was already found 
in 2l|, and through this mechanism ()42p will be obtained also for t < 1. The case of the chain 



defect is analogous. 
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